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ABSTRACT 

Axisymmctric models of the Milky Way exhibit strong interrelations between the 
Galactic constants [the Sun's distance to the Galactic centre (i?o), and the local rota- 
tion speed (0o)], the local stellar columndensity (£*(i?o)) and the shortest-to-longest 
axis ratio of the dark matter halo (q). In this paper we present simple analytical 
approximations that allow for an efficient search through the vastness of parameter 
space, and apply this approximation to investigate the consequences of the uncertain 
gaseous velocity dispersion (a g ) on the constraints imposed by the thickness of the 
Milky Way's gas layer. The extra degree of freedom does not significantly alter the 
conclusions drawn in a previous paper on the shape of the Milky Way's dark matter 
halo. A significant contribution to the total gas pressure by cosmic rays and magnetic 
fields beyond the optical disk is thus ruled out. We find that the Milky Way's dark 
halo is close to spherical if Rq J> 7.1 kpc, while a significantly flattened dark matter 
halo is only possible if our distance to the Galactic centre is smaller than ~ 6.8 kpc. 

Thus, if Rq is larger than ~7 kpc, or 6o <; 170 kms -1 , we can rule out two 
dark matter candidates that require a highly flattened dark matter halo: 1) decaying 
massive neutrinos; and 2) a disk of cold molecular hydrogen. 

It is only possible to construct a self-consistent axisymmetric model of the Galaxy 
based on the IAU-recommended values for the Galactic constants (Rq = 8.5 kpc, 
0o = 220 kms -1 ) in the unlikely case that the effective gaseous velocity disper- 
sion is ~19% larger than observed, and if the local stellar columndensity is less than 
about 18 M^pc -2 . If we assume th at the halo is oblate and a value of S* of 35 ±5 



Mq pc (Kuijken & Gilmore 1989b), we can rule out Galactic models with Rq > 8.0 



kpc and O J> 200 kms -1 . 

Combining the best kinematical and star-count estimates of £*, we conclude that 
£* probably lies between 25 and 45 Mq pc -2 . We find that Kuijken & Gilmore's (1991) 
determination of the columndensity of matter within 1.1 kpc of the plane is robust 
and valid over a wide range of Galactic constants. 

Our mass models show that, largely due to the uncertainty in the Galactic light 
distribution, the dark matter density in the Galactic centre is uncertain by up to three 
orders of magnitude. In the Solar neighbourhood this uncertainty is much reduced: 
our models imply a dark matter density of some 0.42 GeV/c 2 per cubic centimetre, 
or (11 ±5) mM pc" 3 - roughly 15% of the total mass density. 

Key words: Galaxy: structure - Galaxy: kinematics and dynamics - Galaxy: solar 
neighborhood - Galaxy: fundamental parameters - Galaxy: stellar content 
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signatures of the Milky Way's dark matter distribution like 



1 INTRODUCTION 



The observational fact that rotation curves of external galax- 



ies are rather flat in their outer parts ( RuhmFord&^^ 



tano &| vati Gorkom 



aard 198C|; |Bosma 198l|; |Kent 1987|; [Begeman 19871; Paser- 



Broeilri — W?T2r Persic, Salucci 



I- Stcl^ QOfl indicales Ihe presence uf unseen mallei in 
those galaxies^]. A problem which is specific for the Milky 
Way is that we do not know the shape of the Galactic ro- 
tation curve (RC). The slope of the RC depends on the as- 



sumed value s of the Galactic constants (Oiling & Merrifield 



1998b, 1998c). But dark matter is required whatever the val- 



ues of the Galactic constants. A fit to the observed rotation 
curve can be used to yield the stellar mass-to-light ratio Td, 
and the dark halo parameters, albeit with large uncertain- 
ties. Rather than using this fitting procedure, we adopt an 
analytical approximation in which there is only one free pa- 
rameter: the degree to which the disk is maximal (7), from 
which Td and the dark halo parameters follow (see Oiling 
1995 for details). The parameter 7 is preferred for dynamical 
modeling purposes as it is bound between and 1 for a zero 
mass and a full-fledged disk, respective ly. For example, in the 



popu lar "maximum-disk" hypothesis ( van Albada & Sancisi 



1986), the amplitude of the stellar rotati on curve equals (85 
±10)% of the observed rotation speed ([Backett 1997 ), so 
that 7 = 0.85. In contrast, Bottema (1993) used stellar ve- 
locity dispersion measurements to "weigh" stellar disks, and 
concluded that they are sub-maximal, with 7 = 0.63 ± 0.1. 
The situation is similarly indeterminate for the Milky Way: 
the Kuijken & Gilmore (1989, hereafter referred to as KG89) 
model implies 7 ~ 0.5, while more recent models with lower 
rotation speed and shorter scale-lengths can be close to max- 
imal (7 = 0.85 ± 0.1; Sackett 1997). As a result of the un- 
certain stellar mass-to-light ratio, the dark halo parameters 



are very ill determined for most galaxie s ( van Albada 



k Sancisi 1986 



Lake fc Feinswog 1989} piling 1995| ). The 



Milky Way is no different: the combined uncertainty in the 
local disk mass and the stellar scale-length introduces an 
uncertainty of over three orders of magnitude in the central 
dark matter (DM) density of the Milky Way, and about an 
order of magnitude uncertainty in its core radius (§|^;Dehnen 
& Binney 1998). In the Solar neighborhood the situation is 
less dramatic, although no consensus exists on the local vol- 
ume density of dark matter [pz>M or the Oort limit; see §^ 
and Creze el al. (1998) for a recent review]. However, the 
local dark matter density is important for many astrophysi- 
cal problems. For example, if the dark matter comprises ele- 
mentary particles like neutralinos, axions, neutrinos, graviti- 
nos etc., their expected detection rate is proportional to the 
DM density. Likewise, if the dark halo is made up of massive 
compact halo objects (MACHOs), the event rate for gravita- 
tional lensing depends on the integrated dark matter density 
along the line of sight towards the lens. Thus, observational 



micro-lensing ti me scales and optical depths (Gates, Gyuk 
Turner 1995|), and expected neut ralino annihilation rate 



( |Bcrgstrom, Ullio fc Buckley 1998 ) depend on the Galac- 



tic dark matter density distribution and hence the assumed 

values for the Galactic co nstants. 

In a prev ious paper ( piling fc Merrifield 2000, hence- 



forth Paper 1) we determined the dark matter density in 



the Solar neighbourhood at and around R ~ 2Ro to 

infer the minor to major axial ratio (g = c/a), or shape, of 
the dark matter halo of the Milky Way. In the present paper 
we investigate the reliability of several of the assumptions 
made in Paper I, and find that relaxing these assumptions 
does not greatly change the conclusion of Paper I: the shape 
of the dark matter halo of the Milky Way is probably rather 
round. Before going into more detail, let us review some of 
the difficulties which arise when one tries to determine pdm- 
First, large values of the local Galactic rotation speed 
(Oo) result in large DM densities, while low rotation speeds 
require small DM densities. Second, since the shape of the 
Galactic rotation curve depe nds on the value of our dis - 
tance to the Galactic center ( piling fc Merrifield 1998c[ ), 



the amount of dark matter increases with Ro at constant 
Go. And finally, mo re highly flatt ened DM halos have larger 
midplane densities ( Oiling 1995| ). 

We use two sets of observations to constrain the mid- 
plane dark matter density. First, stellar kinematical data 
provides a measure of the total columndensity within 1.1 
kpc of the plane (E^; cf. Kuijken & Gilmore 1991). The 
dark matter density follows after subtracting the luminous 
components and dividing by the scale-height. This method 
yield ambiguous results because uncertainties in the surface 
density of stellar matter (£*) translate into a similar un- 
certainty of the DM density. Thus, low values of S* require 
more dark matter and hence a more highly flattened dark 
halos at constant Ro and Oo. 

Second, the rate at which the thickness of the gas layer 
increases with radius ("flaring") is a measure of pdm- As- 
suming, for now, a hydrostatic balance between internal 
pressure and gravity, it follows that an increasing gas layer 
width is evidence for a decreasing midplane density. Large 
dark matter densities result in thin gas layers, while low den- 
sities yield a thicker gas disk. A larger DM density, due to 
either a larger Oo and/or Rq or a smaller q, results in a thin- 
ner gas layer. However, such a thinner gas layer would also 
occur if the actual gas pressure -or equivalently, the gaseous 
velocity dispersion, a g — is smaller than assumed. Hence the 
significant correlations between the assumed values of the 
Galactic constants, E* and a g and the inferred shape of the 
Milky Way's dark matter halo referred to above. 

In practice, the stellar kinematical data impose corre- 
lations between Oo and q at the Solar circle, while the ob- 
served H I flaring does so at R ~ (2 ± 0.25)i?o • At these 
large radii the stellar disk has vanished so that the poten- 



t See McGaugh & de Blok (1998) for a review of the alternative § Throughout this paper we use cylindrical coordinates with R 

hypothesis that the law of gravity has to be modified instead. the Galactocentric distance, and z the distance of the plane. 
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tial is dominated by the gaseous self-gravity and the dark 
matter. 

Unfortunately, the Galactic constants are ill-determined 



1993 



Oiling fc Merrifield 



I Lvnden-Bell 1986| ; [Reid 
200C), and this has consequences for analyses that 



(Kerr . 
1998a, 

depend on Ro and Go. In an ideal world, parameters such 
as the Galactic constants, the rotation curve [Q(R)], the 
scale-length of the optical disk (/id), the local stellar column 
density, the gaseous velocity dispersion, and so forth, are not 
only measured, but also have normal errors. In that case one 
could determine quantities that depend on these parameters, 
such as the DM density or the halo's shape, by comparing 
the model and observed parameters in a \ 2 sense. Unfor- 
tunately, the quoted errors on Ro and Oo are not normal, 
and the values themselves are not averages in the statistical 
sense, but are rather consensus values with consensus er- 
rors. These arguments lead to the conclusion that the values 
of the Galactic constants and their errors can not be used 
in analyses like maximum-likelihood estimates of a prop- 
erty X which depends on the Galactic constants. Instead 
we urge researchers to investigate the dependence of their 
results on the assumed values of the Galactic constants and 
present their conclusions as functions of Ro and Oo. Given 
the tremendous increase in computing power, such an ap- 
proach is currently more feasible than in the past. Practicing 
what we preach, we follow this ap proach in previous papers 



(Oiling fc Merrifield 1998c, 200C) as well as in the current 



article. 

In Paper I we used both the constraints set by the flar- 
ing of the H I layer as well as the boundary conditions im- 
posed by the local stellar kinematics to infer the halo's flat- 
tening. The procedure works as follows: 1) pick values for 
i?o and Oo and determine the corresponding rotation curve, 
2) create mass models with as ingredients the stellar bulge 
and disk, the interstellar medium (ISM) and a dark halo of 
varying degrees of flattening, 3) pick values for the observa- 
tionally ill-determined scale-length and mass of the stellar 
disk, 4) select a valu e for a g and calculate model flaring 
curves ( |011ing 1995 ), 5) select the model with flattening 
j such that the observed H I flaring at R ~ 2i?o is 
reproduced [this flattening will be independent of the dis- 
tribution of the stellar mass], 6) now vary the stellar mass 
(E*) at Ro , 7) keep track of the dark matter column density 
E^' 1 that needs to be added to E, to match Ej D t, 8) work 
out which halo flattening qx.i is required to generate E^' 1 , 
and finally, 9) select the model that identical gjj j and gj.i 
values. This method is equivalent to finding the zero-point 
of the [qu j (E«) — gi.i(E»)] function. Graphically, this pro- 
cess can be represented as determining the intersection of 
the qui (E») and gi.i(E*) curves (see Paper I, figure 3). 
Note that the qi.i(E*) relation hardly depends on the disk 
scale- lenght. 

In this manner we construct self-consistent mass mod- 
els, with specific values for /id, E* and q, for the selected 
combination of the Galactic constants. Possibly the weak- 
est link in this procedure is the assumed value for o g . In the 
present paper we will remedy this shortcoming of Paper I by 
repeating the procedure outlined above for many values of 



the gaseous velocity dispersion. Picking a value for a g in the 
outer Galaxy is related to the question as to the relevance 
of non-thermal pressure support (due to magnetic fields and 
cosmic ray energy density) of the gas layer. We investigate 
these issues in detail in section ^. 

The surface density of stars in the Solar neighbour- 
hood provides a significant constraint on the determination 
of the halo's sh ape. We therefore review recent determina- 



tions in section 3.1 



Kuijken & Gilmore (1991, henceforth 
KG91) claim that the columndensity within 1.1 kpc of the 
plane is much better determined than the values of the in- 
dividual components. However, their Galactic constants lie 
at the high end of the range we use. Since T,j \ provides 
us with such an important constraint, we consider it pru- 
dent to check KG91's assertion. We indeed confirm KG91's 
finding that Ej D t is relatively well determined over a large 
range in Galactic constants, E* and /id (Appendix |a]). In 
the following section we describe our mass models in more 
detail and determine the density of dark matter in the Solar 
neighbourhood. We summarize and conclude in §ra. 



2 MASS MODELS 

In order to interpret the available data we build axisym- 
metric Milky Way mass models for which we determine the 
model values of the total disk mass, the H I flaring, and 
other parameters. It is well-known that the Milky Way de- 
viates from azimuthal symmetry. However, as we have shown 
in Paper I, and will do so again below, the current observa- 
tional constraints are only barely good enough to rule-out 
the most extreme combinations of Galactic constants. Thus, 
it would be unrealistic to try to incorporate fine-structure 
in the stellar mass distribution due to a bar and/or spi- 
ral structure, especially since these features are not terribly 
well determined themselves. This situation may dramati- 
cally change with the advent of future astrometric satellites 
(e.g., DIVA, FAME, GAIA, SIM) that could determine the 
stellar column density over large parts of the Galactic disk. 

A comparison between the observations and the mod- 
els yields the size, mass and shape of the Milky Way. Our 
models include: a stellar bulge and disk, a gas layer, and 
a dark halo. Some of the relevant parameters of our models 
are tabulated in Table l| We consider models with a range of 
Galactic constants, disk exponential scale-length, stellar disk 
mass, total columndensity, and halo flattening. To calculate 
the exact vertical force law (K 2 (R, z)) at every point (R, z) 
in the Galaxy we integrate over the full mass distribution: 
K„(R, z) = G f °° rdrp(r, 0) dwp(r, w) ± ^L- { with 

s = {r, ui},S = {R, z}, p the total m ass density at (R,z), 
and G Newton's constant of gravity ( Oiling 1995 ). Since 



the calculation of K z and the model flaring curve is rather 
expensive^] we perform these detailed calculations only for 
the limited subset of models listed in Table ^| We determine 



^ The calculation of model gas layer widths for a given combina- 
tion of Ro, ©0) /*d and bulge and disk mass-to-light ratios takes 
about 1.3 hours per (/-values, on a SPARC-10 processor. 
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Table 1. Some parameters used in the model calculations. 



bulge: 



Lib 

h b 
(b 

Pb(0) 

L b ,K,tot 



3.53 
0.61 
667 
10 
15.24 
1.5 10 10 



L &,K P c ~ 

pc 
pc 

L Q,K P c ~ 



disk: 
z e 
h d 
Ld,K(0) 
£*(%) 
P* 

Ld,K,tot 



300 pc 



2, 2.5, 3 



kpc 



2200, 1408, 978 
25-45 
26.5 - 47.7 
5.5 10 10 



L Q,K P c 
L pc- 2 
10" 3 M, 
Lq,k 



(h d =2,2.5,3) 



©PC 



ISM at Ro: 








E H I +H II 


9.25 


M Q pc" 2 






1.80 


M Q pc" 2 




FWHM H j 


410 ±30 


(Ro/7.1) pc 




FWHM H2 


141 ±20 


(-Ro/7.1) pc 




I +H II 


21 ±1.5 


(7.1/flo) lO- 3 


Mq pc" 3 




12 ±1.5 


(7.1/fl ) lO" 3 


M pc" 3 


PISM 


43 ±3.0 


(7.1/-Ro) lO" 3 


Mq pc -3 



the model flaring curves for a much larger range in Ro and 
Go values by applying an approximation (Oiling 1995, Ap- 
pendix D) which we cal ibrate using the models presented in 
Table [2], see section |3.3| for further details. 

2.1 The Mass Components 

2.1.1 Stellar Components 

We base our model for the bulge on Kent's (1992) K-band 
luminosity distribution. The bulge is a modified spheroid 
with "boxy" appearance and density: p b {s) = db Ko{s/hb), 
with s 4 = R 4 + (z I q b ) 4 + Cb ■ The bulge is flattened with axial 
ratio qb- To avoid the singularity of the Tfo-Bessel function 
at s = 0, we include a softening We also truncate the 
bulge exponentially beyond 3 kpc. 

For the disk, we use the standard form of a radially 
exponential disk with central K-band surface brightness 
id,K(0) and scale-length hd- We will consider models with 
value s for the scale- length that are considered to be reason- 
able (Sackett 1997). Furthermore we use the observational 
fact that the total luminosity is better determined than ei- 



ther the central surfa ce brightness or the scale- length (Kent 
Dame & Fazio 1991). Thus, we scale Ld,K(0) such that the 



total luminosity is conserved for each choice of hd- However, 
the numerical values of Ld,K(0) and hd depend on our choice 
of Ro- From Freudenreich (1996, 1998) we find that these 
quantities typically increase with Ro, by 8%x (Ro/kpc - 8). 
Such variations induce a change in total luminosity which 
is a factor of three larger. However, because these correc- 
tions are smaller than the changes resulting from the un- 
certainty in fad, we do not scale Ld,K,tot with Ro- For the 



disk's vertical distribution we use a secant-hyperbolic func- 
tion [sech(2^-)], which is a compromise between the o ften 
used exponentia l and secant-hyperbolic-squared forms (van 



dcr Kruit 1988). For simplicity, the exponential scale-height 
(z a ) of the disk is taken to be constant at 300 pc, intermedi- 
ate between the values suggested by Kent (1992) and Reid & 
Majewski (1993). But note that our results depend neither 
on z c nor on the particularities of the vertical density dis- 
tribution. We use disk mass-to-light ratios (Td,x) such that 
Kuijkcn & Gilmore's (1989b) 2-a range for the local stellar 
columndensity i s spanned. The stellar disk is truncated at 



1996, 1998) 



(Ro + 4.5) k pc ( [Robin, Creze fc Mohan 1992 



Freudenreich 



2.1.2 Interstellar Medium Components 

We now turn to the contribution to the mass models from 
the interstellar medium. Our location inside the Milky Way 
means that distances to diffuse components like the ISM 
are based on a kinematical model. Thus, important proper- 
ties like the full width at half maximum (FWHM), volume 
density, and total mass of the ISM depend on the Galac- 
tic constants and the Milky Way's rotation curve. Thus, we 
re-determine the atomic and molecular gas distributions for 
each choice of Ro and ©0 . While the columndensities at 
fractional radius R/Ro are independent^ of the choice of 
Ro & Oo, the number density, the thickness and the to- 
tal gas mass of the Galaxy are not. Observationally, we 
can currently do no better than determining the gaseous 
column-density (E 9 ) at fractional radius R/Ro- Likewise, 
since the thickness measurement returns an angular size, we 
can only determine FW HM/_Rp at fractional radius R/Ro 
[FWHM/flo = C,R/R , ( |Binney fc Merrifield 1998| ), their 
figure 9.25]. Physical length-scales for the thickness and 
Galactocentric radius can be assigned after choosing a value 
for Ro. 

For the inner Galaxy, the H I colum ndensity was d eter- 
mined from the midplane volume den sity ( Burton 198£ ) and 
the observed thickness of the layer (Malhotra 1995). The 
H I columndensities for R > Ro were taken from Wouter- 
loot et al. (1990). The H2 columndensities for the inner and 
outer Galaxy were copied from Bronfman et al. (1988) and 
Wouterloot et al. (1990), respectively^*]. For the radial pro- 



I From Binney & Merrifield's (1998) equation 9.12 it follows that 
the distance-dependent function that appears in the exponential 
term only depends on r = R/Ro, so that one can rewrite the 
observed brightness temperature at line-of-sight velocity M; os as 
Ro J drF(r, ...)n(r)H(ui og ,r, ...) with F a function that results 
from the transformation from line-of-sight distance to r, and H 
the velocity-dependent function. If H varies more quickly with 
r than the number density n(r), then the product Ro X n(r) is 
approximately constant. Thus, as Ro is increased, the volume den- 
sity goes down. Because the gas-layer width increases with Ro, 
the vertical column-density at scaled distance R/Ro is indepen- 
dent of the Galactic constants. See also Bronfman et al. (1988). 
** A graphical representation of the derived gaseous surface den- 
sity distributions can be found in a related paper on the Galactic 
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files of the ISM we neglect the columndensity due to the 
other phases of the ISM. However, in order to properly take 
into account all baryonic contributions to the local c olumn- 
density , we include 1.4 M pc -2 of ionized hydrogen (Kulka- 
rni & Heiles 1987 ) at the Sol ar position only. We furt her 
include 23.8% helium by mass ( Olive fc Steigman 1995 ) to 
compute the surface densities. 

As in our previous papers, we do not distinguish be- 
tween the various phases of the atomic Hydrogen but rather 
assume that the warm neutral medium has a kinetic temper- 
ature equivalent to the "temperature" associated with the 
bulk motions of the clumped, cold neutral medium. Note 
that the cold medium is likely to be absent beyond the stel- 
lar disk ([Paper l|; |Braun 199% 



2.1.3 The Dark Matter Component 



As is commonly done (e.g., van Albada & Sancisi 19861 


Kent 1987; 


Bcgcman 1989; Lake & Fcinswog 198E ; Sackett 


& Sparke 199C ; 


Brocils 1992 


; Oiling 1995; Oiling 1996b), 



we model the dark matter density distribution as a non- 
singular isothermal spheroid with flattening q, and a density 
distribution given by: 



Ph(R,z;q) = p (q) 



R 2 c (q) 



R 2 c (q)+R2 + (z/q)< 



(1) 



where the halo's core radius (R c ) and central density (po) 
depend on the flattening in such a way that the family of 
density distributions ph(q) have a rotat ion curve that is es- 
sentially independent of q ( Oiling 1995 ). Figure [l] shows the 
radial distribution of dark matter for three values of Ro, for 
the range in model parameters as listed in Table [| Each 
choice for hd and E, results in a different distribution. The 
largest DM densities are obtained for large /id's and small 
E„'s. Even though the central DM density is uncertain by 
over three orders of magnitude, the local dark matter density 
is determined rather well: we find 



PdmCBo,9o) _ 11.5 + 3.8 x (Ro - 7.8) ± 2 
lCT 3 M pc- 3 ~ <?(26.7/fi ) 2 



(2) 



with f2o = Oo /Ro, and where Ro is in kpc, and Oo in kms -1 . 
The (~ 25/q)% uncertainty arises as a result of the un- 
certainty in Ej t and E*. For example, taking Ro = 7.1 
kpc, q = 0.71, and using the parameters tabulated in Ta- 
ble 2, we find pdm(Ro) = 10.5 mM pc~ 3 (1 mM pc~ 3 = 
10" 3 M pc~ 3 ). Thus, the R = 7.1 model value for the 
DM matter density compares well with the values given by 
uation (0) as well as with the observational determination 



(§3-1) 



Oort constants (Oiling & Merrifield 1998c). The H2 column 



densities have been re-scaled from the original sources assuming 
N(K 2 )/W(CO) = 2.3 10 20 cm" 2 (K kms" 1 )- 1 . We present the 
radial increase of the thickness of the gas layer in Paper I. The 
thickness and columndensity of the H I and H2 are tabulated in 
Appendix O 



CO 

I 

a 

Oh 

© 



o 
o 



R =7.1; q=0.75 
R„=7.B; q=1.00 
R =8.5; q=1.00 




Galactocentric radius [R ] 

Figure 1. The dark matter midplane density distribution calcu- 
lated using equation (flh for three values of Ro an d q. The densities 
scale approximately as 1/q (Oiling 1995) The error bars represent 
the full range for the case corresponding to E* = 35±5 M pc — 2 , 
and /id=2 to 3 kpc. Notice that although the DM density at the 
Solar circle is relatively well determined, the central density is 
uncertain by more than three orders of magnitude. Note that we 
have offset the radial coordinates of the i?o=7.1, and 8.5 kpc 
models from their true positions (i.e., the iJo=7.8 kpc points) to 
avoid overlapping symbols. 



2.2 Other Dark Matter models 

Of course, other mass models can be constructed w hich rep- 



resent the radial dark ma t ter density distribution ( 



Albada K 


z Sancisi 


1986 


Dehnen & 


1 Binney 


1998) 



_ 



Navarro, Frenk & White 1996 



However, all viable mass models 
must share the property that they reproduce the Galactic 
rotation curve. For a round halo the vertical force approx- 
imately equals z/R times the radial force. Since the radial 
force is approximately the same for all models which repro- 
duce the observed rotation curve, the ensemble of possible 
models also have approximately equal vertical forces , inde- 
pend ent of the exact radial DM density distribution ( Oiling 
1995). In a flattened halo with the same rotation curve, the 



DM densities are roughly proportional to 1/q, independent 
of the radial mass distribution [cf. eqn. (Q)]. Thus, our analy- 
sis will not be seriously compromised by restricting ourselves 
to one particular DM density distribution. 
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Table 2. Disk & halo parameters. A log for a representative sample of models for which exact flaring curves are calculated. The meaning 
of the columns is as follows: 1) 7, the degree to which the stellar disk is maximal; 2) scale-length of stellar disk (kpc); 3&4) bulge & disk 
K-band mass-to-light ratios; 5&6) stellar columndcnsity & total columndensity within 1.1 kpc of the plane (Mq pc — 2 ); 7&8) core radius 
(kpc) and central density (uiMq pc — 3 ) of the dark halo, assuming q=l (multiply by s=s 1/q to scale to different values of q, for a more 
exact scaling, see Oiling 1995); 9&10) dark halo flattening as inferred from the H I flaring (gjj j ) and the local columndcnsitie s (gi.i). 
For a given set of (Ro,@o,h^), the gjj j errors for the Milky Way are similar to those for extra-galactic systems (Oiling 1996a) . 
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yiX.l 

^tot 


Rc po(q = 1) 


m 1 


91.1 


i?o=6.8, O =182 



0.80 2.0 

0.86 2.0 

0.92 2.0 

0.64 2.5 

0.70 2.5 

0.75 2.5 

0.56 3.0 

0.60 3.0 

0.65 3.0 



i?o=7.1, e =189 



0.82 2.0 

0.89 2.0 

0.95 2.0 

0.65 2.5 

0.71 2.5 

0.76 2.5 

0.56 3.0 

0.61 3.0 

0.63 3.0 



R =7.8, e =207 



0.92 2.0 

0.72 2.5 

0.78 2.5 

0.84 2.5 

0.59 3.0 

0.61 3.0 

0.68 3.0 



i?o=8.5, ©o=227 



0.76 


2.5 


0.60 


0.670 


30.66 


78.68 


6.81 


39.1 


2.922 ±0.863 


1.169 ±0.252 


0.81 


2.5 


0.60 


0.761 


34.83 


80.43 


8.34 


28.9 


2.512 ±0.608 


1.228 ±0.279 


0.86 


2.5 


0.60 


0.857 


39.26 


81.28 


11.10 


19.9 


1.931 ±0.366 


1.256 ±0.306 


0.59 


3.0 


0.85 


0.486 


27.23 


76.53 


5.09 


60.2 


1.424 ±0.156 


1.165 ±0.193 


0.67 


3.0 


0.85 


0.627 


35.11 


82.49 


6.20 


43.2 


1.498 ±0.183 


1.364 ±0.213 


0.75 


3.0 


0.85 


0.785 


44.00 


88.53 


8.07 


28.9 


1.475 ±0.197 


1.641 ±0.159 



0.55 


0.422 


30.21 


62.48 


2.39 


73.8 


0.630 


±0.107 


0.561 


±0.175 


0.55 


0.488 


34.92 


66.02 


3.54 


35.6 


0.687 


±0.113 


0.664 


±0.226 


0.55 


0.558 


39.96 


68.72 


6.09 


14.6 


0.699 


±0.120 


0.753 


±0.295 


0.55 


0.327 


29.58 


63.56 


1.32 


246.3 


0.687 


±0.152 


0.620 


±0.183 


0.55 


0.392 


35.39 


68.84 


1.86 


125.1 


0.717 


±0.134 


0.826 


±0.268 


0.55 


0.449 


40.63 


73.40 


2.45 


73.0 


0.674 


±0.127 


1.061 


±0.371 


0.55 


0.306 


30.21 


64.34 


0.70 


853.9 


0.740 


±0.181 


0.651 


±0.194 


0.55 


0.351 


34.68 


68.6 


1.00 


421.0 


0.770 


±0.159 


0.820 


±0.261 


0.55 


0.412 


40.70 


74.18 


1.44 


201.9 


0.762 


±0.137 


1.114 


±0.382 



0.45 


0.482 


29.71 


61.94 


2.73 


62.7 


0.625 


±0.132 


0.544 


±0.165 


0.45 


0.568 


35.00 


65.66 


4.53 


25.5 


0.624 


±0.116 


0.636 


±0.214 


0.45 


0.647 


39.87 


66.74 


9.77 


8.6 


0.679 


±0.104 


0.603 


±0.239 


0.55 


0.370 


29.66 


62.60 


1.13 


340.1 


0.727 


±0.159 


0.581 


±0.173 


0.55 


0.441 


35.39 


67.76 


1.66 


158.2 


0.738 


±0.161 


0.765 


±0.250 


0.55 


0.506 


40.55 


72.26 


2.25 


86.4 


0.823 


±0.223 


0.999 


±0.356 


0.55 


0.333 


29.73 


61.82 


0.20 


9884.5 


0.735 


±0.181 


0.550 


±0.168 


0.55 


0.395 


35.28 


67.10 


0.55 


1343.6 


0.703 


±0.158 


0.734 


±0.243 


0.55 


0.421 


37.63 


69.32 


0.71 


802.8 


0.786 


±0.210 


0.843 


±0.290 



0.45 


0.687 


29.80 


64.60 


9.51 


15.4 


1.003 ±0.119 


0.616 ±0.162 


0.75 


0.492 


29.80 


64.76 


4.59 


47.9 


1.744 ±0.553 


0.963 ±0.240 


0.75 


0.577 


34.97 


71.00 


5.81 


32.6 


1.616 ±0.443 


1.065 ±0.260 


0.75 


0.669 


40.56 


77.66 


8.01 


20.6 


1.461 ±0.337 


1.191 ±0.331 


0.75 


0.430 


30.46 


66.54 


2.84 


103.4 


1.623 ±0.451 


0.973 ±0.244 


0.75 


0.460 


32.56 


68.06 


3.02 


91.7 


2.019 ±0.758 


1.026 ±0.260 


0.75 


0.572 


40.46 


77.54 


3.81 


59.3 


1.810 ±0.753 


1.327 ±0.364 



3 CONSTRAINING THE DARK MATTER 
DENSITY 

In this section we present two simple analytical models to 
illustrate the existing correlations between Ro, Oo, S tot , E«, 
q, and a g outlined in the Introduction. These models can be 
used to show how the local stellar kinematics and the H I 
flaring constrain the dark matter density in the Solar neigh- 
bourhood and at ~ 2 Ro- But first we investigate how well 
the local stellar columndensity is determined and how accu- 



rately this determination constrains the dark matter density 
in the Solar neighbourhood. 

3.1 The local stellar columndensity: a constraint? 

The Milky Way is a unique galaxy in that we can, at least 
in principle, determine the local columndensity of stars di- 
rectly. Once £* is accurately determined, it is possible to 
establish to what degree the Milky Way disk is maximal, 
which would provide an important benchmark for external 
galaxies. Unfortunately this benchmark is not yet available, 
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as the values for E* reported in the literature range from 26 
to 145 M pc -2 . 

Two basic techniques have been employed to determine 
£«. The direct method involves converting star counts as 
a function of Galactic coordinates and magnitude to in- 
situ mass densities. This method is somewhat hampered 
by uncertainties in the conversion from luminosity to mass, 
completeness problems in the Solar neighborhood, and bi- 
nary corrections at large distances. Gould, Bahcall & Flynn 
(1997; hereafter referred to as GBF97) used deep HST star 
counts of M-dwarfs at great heights above the plane in com- 
bination with a local normalization to infer a stellar column- 
density of only 25.8 ±3.8 Mq pc -2 . The second, kinematical, 
method employs the interrelationship between the potential, 
the vertical density distribution, and the variation of the ve- 
locity dispersion with z. Many authors have employed this 
method, yielding a large range in inferred values for E*. For 
example, Bahcall (1984b) found that the inferred value of 
E* depends significantly on the assumed vertical distribu- 
tion of the dark matter: ~40, ~52, ~68, and ~145 Mq pc -2 
for dark matter distributions that resemble the gaseous disk, 
an isothermal halo, the thin stellar disk, or the thick stel- 
lar disk, r espectively. Other authors find values as low as 3 5 



M Q pc" 2 (Kuijken fc Gilmore 1989b ; Flynn fc Fuchs 1994) 



Methods that combine the two pr i mary methods exist as well 
(Bienayme, Robin & Creze 1987; Creze, Robin & Bienaymc 
198S|) . In Appendix |a| we simulate the kinematic determi- 
nation of the stellar columndensity with the aid of Galaxy 
models -taken from Table |^- for which we know the exact 
force law. Our results are in complete agreement with the 
findings of previous authors: the "vertical disk-halo conspir- 
acy" can only be resolved if high-z stellar kinematical data 
are included and/or if additional assumptions are made (e.g., 
a "reasonable" value for the DM de nsity) . Typical v alues for 
t he mass of the stellar disk are 52 (Bahcall 19841 ), 35 ± 5 



(iKujiken fc Gilmore 1989b|) and 37 ±13 Mq pc~ 7 ( |Flynn & 



Fuchs |L994| ). 

Considering the uncertainties in the kinematical esti- 
mates of the total mass of the disk, it might be preferable to 
use the direct star count method to determine E*. The lat- 
est results by GBF97 imply E»=25.8 ±3.8 A/ Q pc -2 . How- 
ever, the local space density of stars found by these authors 
[(33.2 ± 8.6) vciMq pc -3 ] is somewhat lower than the values 

reported in the literature [(43 ±15)^]mAf0 pc~ 3 ]. This sug- 
gests that GBF97 may have under-estimated E, by a factor 
of 1.4 ±0.2, and so we find E* = 36 ±5 Mq pc" 2 : the star 
count method yields a value for £„ that is remarkably close 
to the kinematical estimates. 

On the other hand, the total matter density in the So- 
l ar neighbourhood as inferred from recent Hipparcos data 



(Creze et al. 1998) of 76 ±15 hiA/q pc favours a local 



stellar volume density which is even smaller than reported 
by GBF97. After subtracting the density of the ISM (see 



tt T he average of: p* = 46 ( Widen 
1984 4 fcahcall, Flynn fc Gould 
et al. (1998). 



1974[ ), 50.8 ( |Bahcal] 
199^ ), 45 wlMq pc" 3 [CrTSe 



Table |l|) and the local DM density [eqn. (0)] we find a local 
stellar density of ~ 21 ±15 uiMq pc -3 . Note that even lower 
values for p* result if large values for the Galactic constants 
and/or Slo are chosen, while smallish Galactic constants and 
Qo increase the local stellar density. 

To summarize, a "reasonable" value for E* might be 
35 Mqpc -2 , and we suggest a "consensus" error of 10 
Mq pc -2 . However, since the differences between the various 
E„ estimates are not random but systematic, one should not 
interpret the reasonable E* value and its error in a statis- 
tical sense. That is to say, one can not construct likelihood 
contours based on a reasonable average value and a con- 
sensus error bar. Any attempt to do so would be an over- 
interpretation of the available data. 

A "reasonable" value for the local dark matter column 
can be obtained by subtracting the columndensities of the 
stellar and ISM distributions: E h = (24.5 ±11) Mq pc -2 
within 1.1 kpc of the plane. This dark matter column 
amounts to an average local dark matter density of 11 ±5 
mM pc~ 3 . 

3.2 The connection between dark and luminous 
matter in the Solar neighbourhood 

In the previous section we have seen that the contribution 
of the stellar and dark matter components to the local disk 
mass cannot be clearly segregated. Thus the observed value 
for Et i implies that E* and Eh are highly correlated: a low 
stellar column implies a large amount of dark matter, and 
vice versa. Since the dark matter density depends on the 
amplitude of the rotation curve and the halo's flattening, 
these parameters are in turn related to E*. Below we in- 
vestigate the relations between T,l \, Oo, E, and q in some 
detail. These relations are independent of the H I flaring. 

Integrating equation (|l|) with respect to z and applying 
equation (A4) from Oiling (1995), we find the columndensity 
of dark matter within z of the Galactic plane as a function 
of q: 



mi) 



*j0mi - z/q 



2-kG y/Rl(q) +R 2 arccos < 



■ arctan 



z / ( l 



v/i? 2 (g)±i? 2 



(3) 



. (4) 



with Vh.oo the asymptotic rotation velocity of the round dark 
halo. For each of our self-consistent mass models we can 
calculate E^(g) and find that it obeys a simple power law 
relation: 



£i- 1 (g) + E,±E g 



(5) 

w £&(g = l) x g -». (6) 

The dependence of T*l \ on the other parameters of the 
model parameters are combined into the index p (with 
p~ 0.05 -0.45)0 The halo flattening inferred from the lo- 

tt The values for E^(g), £^(? = 1) and p are different for 
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cal stellar kinematics (gi.i) is then found by equating E 4 „t(g) 
with the observed value and solving equation (^) for q 



= 



71 ±6 



(7) 



Equations (|4|), (JfJ) and (|7|) can be combined to find the de- 
pendence of gi.i on the Galactic constants. Neglecting R c 
and equating 14, oo with Oo, we find: 



9i.i 



cQg/fl 2 , + S, + E s 

yl.l 



VP 



(8) 



where c is a constant. Qualitatively, this is the functional 
dependence we expect: increasing Oo at constant Rq would 
increase the required amount of DM, which has to be coun- 
teracted by increasing gi.i to retain the same Eto t . Likewise, 
increasing Ro would place the Sun in a lower density region 
of the halo, which needs to be compensated by decreasing 
qi.i. Furthermore, given the small values of p, qi.i depends 
very strongly upon the Galactic constants. 

In Appendix |b| we present some other correlations be- 
tween the Galactic constants, E», ©o and qi.i. For example, 
equation (B4) reveals the linear relation between gi.i and E* 
for q > 0.15. Note that models with different scale-lenghts 
follow the same gi.i(E») relation, albeit that different /la's 
yield different gi.i's for a given stellar column density^} 
These examples show that the E* t constraint implies highly 
flattened DM halos for small values of Oo. 

3.3 Constraints from the thickness of the gas layer 

In an idealized picture, the equilibrium thickness of the gas 
in the Milky Way depends only on the gas "temperature" 
and the form of the potential in which it has settled. Thus, 
the observed FWHM of the gas layer can be used to con- 
strain the potential of the Galaxy. In Paper I we present ev- 
idence that the interstellar medium beyond the optical disk 
comprises only a single, iso-thermal component. We there- 
fore adopt the same assumption in the analysis below. In this 
section we will expand our analysis to include the effects of 
non-thermal pressure gradients. 

If the gas layer is in a steady state and we assume that 
only "thermal" motions of the gas contribute to the pressure 
in the ISM, the thickness of the gas layer follows from the 
equation of hydrostatic equilibrium: 



dz 



p s (z) K z (z) 



(9) 



where a g , p g (z), and K z are the gaseous velocity dispersion 

each combination of the Galactic constants, the stellar disk mass 
and scale-length, p can be determined from a fit to equation (JgJ) , 
where we determine Ej ^(g) for several models with varying q's. 
Note that we also use equation (^j) to extrapolate into the prolate 
regime, with the exponent p as determined for the oblate models. 

The reader can verify this by plotting the qi.i values as a 
function of S* from table bl 



and volume density, and the vertical force per unit mass, 
respectively. The vertical force is commonly determined us- 
ing a "local approximation" which K z follows from the the 
Poisson equation and the local mass densities and the ra- 
dial gradient of the rotation curve. However, this approach 
is known to fail in regions where either the mass densities 
or the rotation curve, or both, have steep gradients. Fur- 
ther, the local approach neglects any variation of the circular- 
speed with height above the plane. These problems can be 
overcome by calculating K z from the global mass distribu- 
tion (the "global approach," Oiling 1995). As mentioned in 
section ^, we combine the exactness of the global approach 
with the speed of the local approximation for optimal re- 
sults. 

If we employ the local approximation the dependencies 
between the various parameters of the model become appar- 
ent. For example, when the potential is dominated by mass 
component i, equation ^ can be solved for the thickness of 
the gas layer: 



FWHMi oc 



^/47^Gp~ 



(10) 



where the proportionality constant depends on the vertical 
density distribution of the dominant contributor to the lo- 
cal vertical potential. We copy some relations from Oiling 
(1995). In case the gas is fully self-gravitating, the width is 
given by: 



FWHM e 



0.158^2- 

E 



(11) 



If the potential is dominated by an iso-thermal stellar disk 
with sech 2 scale-height zq (zo = 2z e ), the thickness of a gas 
layer would be: 



FWHM, 



0.51a 



99.2 



£0,0.6 



E, 



(12) 



or 510 parsec at the Solar circle (a g , zo and E* are expressed 
in units of 9.2 kms" 1 , 0.6 kpc and 35 Mq pc -2 ; see also, van 
der Kruit 1988). And if the dark matter halo dominates the 
potential, we find: 



FWHM h 



13.5 q 



1.4 + q 14,, 



Ki + R2 ■ 



(13) 



where R c ,i is the core radius of the equivalent round halo 
( |011ing 1991: ). All distances and widths are in kpc, and 
all velocities in kms" 1 . The thickness of the gas layer in 
the combined potential of several mass components can be 
solved analytically in the form of an integral equation, but 
requires an iterative solution procedure (Oiling 1995, Ap- 
pendix C). However, a further approximation is possible. 
Following Oiling (1995, Appendix D) we use: 



FWHM 



E 



FWHM 



(14) 



where the weighting factors Wi reflect the relative impor- 
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tance of component i in the region where the gas reside sQ 
Solving equation [TiJ for qjj j , taking Vh,ao ~ ©o and ne- 
glecting R c ,i, we find: 



<?H I (2%) 



1.4 



')( 



2Ro 
w h 



(15) 



with £{R) = FWHM obs /i? (~ 0.07 for the outer Milky Way), 
and where S s , w g , and Wh have to be evaluated at R = 2i?o- 
From this equation we infer that the contribution of the 
self-gravity of the gas depends strongly on the gaseous ve- 
locity dispersion: low a g values result in a more negative self- 
gravity contribution to the denominator, and hence leads to 
a larger inferred gjj t . As we have already mentioned, the 
self-gravity of the gas is an important component as it re- 
duces the width of the gas layer by approximately 45%. A 
further simplification can be made by neglecting the gaseous 
self-gravity as well, we find: 



1R I 



9H I 



1.4[C(-R)6 ] 



13.5a,, 2 - [((R)Qo] 2 ' 
9o 2 



(16) 
(17) 



where the second line arises because the first term in the de- 
nominator dominates. The errors in the halo flattening and 
rotation speed are related through equation (CI), where we 
have neglected the contribution of the velocity dispersion 
error. The dependence of gjj j on Go and a g are indeed 
as outlined in the Introduction. Also note that, due to the 
quadratic nature of the proportionalities, the observed flar- 
ing and the small allowed range of qji \ constrain 0o and 
(j g rather tightly. 

Comparing equations (^) and (^) we see that the con- 
straints on the halo shape arising from the local stellar kine- 
matics and the H I flaring the have rather different Ro and 
Oo dependencies. Thus, it is indeed possible to learn more 
about the Galactic dark matter distribution by combining 
these two constraints. Further, unlike the stellar kinematical 
method, the constraints from the H I flaring is independent 
of E* since it arises at R ~ 2Ro , beyond the truncation of 
the stellar disk. 

Equations (^j|)-([l7|) only serve to illustrate the depen- 
dence of 5fj j on the model parameters since several im- 



H Sec Oiling (1995) for details. For the Milky Way we use: 
w g = ((1.24£ fl + 1.55W h x (p h + p r ))/(Zg + W h x (p h + p r )), 
w h = (1.0 + 0.35 x exp (-LW^t/E,)), and w r = (1 + w h )/2. 
ph is the local dark matter density, the "rotational" density, p r (= 
2 ~g V r 1 dV dR t )' ar i ses due to the slope of the rotation curve. 
W h equals FWHM h /2.35. For example, at R~ 2R , S 9 ~ 7.7 
and (E h ,£ r ) = (3.6, -0.30), (5.3, -0.17), (6.0, -0.37) M pc" 2 
for i?o=7.1, 7.8, and 8.5 kpc, respectively. These values lead to 
weights that are almost independent of Ro: w g ~ 1.4, ~ 
1.2, uv ~ 1.1. These weightings embody the calibration of the 
local approximation for the outer Milky Way. 



portant aspects are treated too simplistic. First, neglecting 
R c and equating Vh,oo with Oq is only seldomly warranted. 
Second, depending on the slope of the rotation curve, an 
error of order ±20% is made in the inferred 5hj_. We there- 
fore do not recommend using equations (|l5|)-([l^) "as is" to 
determine qu j directly. In all our calculations we employ 
equation (114) where we take all mass components fully into 
account, without any further simplifications. For many com- 
binations of model parameters, a round halo is too dense to 
explain the observed flaring (cf. Fig. Albeit not entirely 
correct, we will determine the halo's contribution to the po- 
tential using equation (|l^) in those cases. 

4 THE EFFECTS OF A LARGER GASEOUS 
VELOCITY DISPERSION 

Equations (|^)-((l^) show that gjj i is a function of a g , 
Go and Ro- On the other hand, the halo flattening inferred 
from the local stellar kinematics is a function of the Galactic 
constants and £*. In a self-consistent model, qi.i = gjj i , so 
that strong inter-relations are imposed among the currently 
ill-determined values of Ro, Oo, E* and a g . 

Below we investigate how the inferred halo flattening 
and stellar columndensity of a self-consistent model depend 
upon our choice of a g . The value of a g can have significant 
effects on the inferred values of q and E*. For example, in 
Paper I we assumed that the true velocity dispersion equals 
9.2 kms -1 , and found an upper limit to the local rotation 
speed of about 190 kms -1 . Further, models with the IAU 
recommended Galactic constants have prolate dark halos 
(q ~ 1.9) and require a rather high local stellar columnden- 
sity of ~55 Mqpc -2 . Increasing a g would bring q and E* 
down to more acceptable levels, and it might thus be worth- 
while to treat a g as a free rather than as a fixed parameter. 

However, observations of external galaxies imply that 
the velocity dispersion declines slightly in the radial range 
over which the Milk y Way's flaring has been measured, by 
a factor 1.12 ±0.12 (jvan dor Kruit fc Shostak 1984 |Dickey,' 



1996b 



Hanso n fc Helou 199C| ; |Kamplmis I' >'>■'.[ |C6tc 1995| ; piling 
Sicking 1997[ ). rmore, there is some evidence 



that a g in the outer Galaxy equals the value inside the Solar 



circle (Blitz & Spergel 1991). Also, one would expect that 



a change in gaseous velocity dispersion would be reflected 
in a change of the residual motions of young stars with re- 
spect to the mean streaming field. Such changes have not 
been observ e d, neither in B stars nor in Cepheids ( Brand & 
Blitz 1993| ; front, Queloz, Bratchi fc Mayor 199?| ). Thus, 



an increase in gaseous velocity dispersion beyond the Solar 
circle is not likely. 

4.1 Non-thermal pressure terms 

In Paper I we argued that the ISM beyond the optical disk 
comprises a single, iso-thermal component (the warm neu- 
tral medium). We now investigate the possibility that non- 
thermal pressures have to be included in the hydrostatic 
balance. 

It is estimated that in the Solar neighbourhood thermal 
motions, cosmic rays, and magnetic fields contribute about 
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equally to the interstellar pressure ( |5pitzer 1976 ; Kulka- 
rni & t eiles 1987). However, note that the gas-layer width 

not just the pressure 



depends on the pressure gradient, not just the pressure. In- 
corporating the non-thermal pressure terms, the equation of 
hydrostatic equilibrium can be written as: 

/ \ K- i \ dPa , dF>B , dPc dPtot 

Pe (z)K z (z) = — + — + — = — (18) 



the thermal gas pressure, Pb — B /8iv 



with Pg = <J g 2 p g (z) 
the magnetic field pressure, and Pc = l/3Uc is the pres- 
sure due to cosmic rays with energy density Uc- Taking the 
scale-heights of the kinetic, magnetic and cosmic ray energy 
density to be zg,zb, and zc, we write the pressure gradients 
as ^ = f£, = and 



dP C 

dz 



z c ' 



4.2 Gas-Layer Support at Ro 

In the Solar neighborhood the scale- heights zb and zc 
are Kb and kq time s larger than the gaseous scal e -height 



1KB 



10. 



4; (s.g 



Kulkarni & Heiles 1987 



Rupen 



1991)1- Following Spitzer (1978) we assume that the non- 



thermal pressure terms are proportional to the kinetic pres- 
sure: Pb = o-bPg and Pc = acPc, with as ~ 0.25 and 
etc ~ 0.4 in the Solar neighborhood. Using these parame- 
terizations we find: 



p e (z) K*(z) 



°9 2 ( 



1 + 



— I a ° \ d P° 



KB 



kc J dz 



(19) 
(20) 



where we have lumped all terms contributing to the hydro- 
static balance into a single unknown, the effective veloc- 
ity dispersion a'„. With the above simplifications and equa- 
tions (y_0|) and (|16[), it becomes possible to estimate the ef- 
fects of non-thermal pressure support on the thickness of the 
gas layer and the inferred shape of the dark matter halo: 



FWHMgbc 



FWHMg 



Igbc 



^ otB ate 

KB KC 



1.06 



'la 



0.89. 



(21) 
(22) 



where the subscripts GBC indicates that the gaseous, mag- 
netic and cosmic ray terms are taken into account in the 
hydrostatic balance of the gas layer. Thus, neglecting the 
non-thermal contribution to the pressure balance in the lo- 
cal ISM leads to a gas layer width that is under-estimated by 
a few percent only. Similarly, a slightly more flattened halo 
is required if non-thermal pressure support were important 
in the Solar neighborhood. 



4-2.1 Non-Thermal Support, or Not? 

The above statements are at odds with the current paradigm 
[cf. Binney & Merrifield (1998), problem 9.7] which states 
that there exists significant non-thermal pressure support 



of the gas layer in the Solar neighborhood. The argument 
in support of this paradigm is that models without non- 
thermal pressure support under-predict the observed gas- 
layer widths. However, the model predictions depend sensi- 
tively on the values of the Galactic constants. Our models 
with small values for Ro and Oo show only a small discrep- 
ancy between the model and observed widths, while a large 
width difference exists for models with IAU-standard Galac- 
tic constants (Paper I, figure 5). Furthermore, the gas- layer 
width depends sensitively on the local stellar column density 
as well as the functional form of the vertical density distri- 
bution. For example, the gas layer is almost 40% thinner in 
case the stellar ve rtical density distrib ution is exponential 
rather than sech 2 (van der Kruit 1986). 



In fact, if the potential in the Solar neighborhood were 
fully determined by the stellar disk, and the H I were truly 
iso-thermal, then the predicted width via equation ^| over- 
predicts the observed width (cf . table |l|) . The other contrib- 
utors to the potential decrease the model width slightly (cf. 
cqn. |l4|), in better agreement with the observed width. We 
conclude that the paradigm of significant non-thermal pres- 
sure support in the Solar neighborhood is inaccurate and 
that the observed H I width at the Solar circle is consistent 
with "thermal" pressure support only (cf . equation ^l|) . 

4.3 Gas-Layer Support at 2Ro 

Since we only employ the H I thickness measurements at 
large radii as a constraint in our mass models, let us try to 
guess the value of o' g at ~ 2Rq. 

The gas density in the outer Galaxy is about four times 
smaller than at Ro as a result of the decrease in column- 
density (factor 2) and the increase in thickness (factor 2). 
To make a conservative estimate as to the importance of 
the non-thermal terms at large distances, we assume that 
the magnetic field strength, the cosmic ray energy density 
and their vertical gradients equal the values in the Solar 
neighbourhood. Furthermore, we assume that o g remains 
unchanged. With these assumptions, the a values at 2Ro are 
four times larger than at R = Ro, while the k values are 
decreased by a factor two. Thus, the effective velocity dis- 
persion and the inferred halo shape are strongly affected: 



4a s (i? ) . 4a c (-Ro 



+ ■ 



qcBC 



'la 



kb(Ro)/2 k c {Ro)/2 



^ -0.5. 



V2 (23) 



(24) 



On the basis of this over-estimation, it appears that the 
gas layer is significantly out of thermal equilibrium at large 
radii. Do we expect this to be the case? Two effect reduce 
the importance of the non-thermal contribution. First, the 
B-field is probably co-spatial, "frozen in", with the ionized 
part of the ISM which' scale-height has most likely doubled 
at R — 2Ro (like the HI). Thus, the Solar neighbourhood 
value for kb(2Ro) should be used in equation ( J23|) . Second, 
radio-continuum measurements of external galaxies suggest 
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that high cosmic ray fl uxes are closely associ ated with the 
sites of star formation (Bicay & Helou 1990). Since active 
star-formation ceases to exist beyond ~ 1.5i?o, the Qc-term 
in eqn. vanishes. Thus, in a more realistic treatment of 
the non-thermal pressure terms, their effect on the vertical 
equilibrium of gas at large radii is much reduced: 



qoBC 



qo 



2E\, 



2R„ 



ia B {Ro) 



1 + 



k b (Rq) 



^ -0.9. 



TT ~ 1.05 (25) 
(26) 



Further, if the magnetic field strength decreases with Galac- 
tocentric radius, the above estimates should be even closer 
to unity. Thus realistic estimates of the importance of non- 
thermal pressure terms in the equation of hydrostatic equi- 
librium indicate that the the thickness of the H I layer is not 
much affected, neither in the inner, nor in the outer Galaxy. 

And finally it is worth mentioning that there exist both 
theoretical and empirical evidence that non-thermal pres- 
sure support is small even when simplified calculations in- 
dicate that they may dominate dPc/dz. From a theoretical 
perspective one finds that the vertical balance is only af- 
fected if the magnetic field is horizontal ly stratified, a con- 
figuration that is unstable (Parker 1966), and which is thus 
not expected to exist on a global scale in the Milky Way. In 
fact, 3D simulations of the growth of the Parker instability 
show that the gas is almost entirely supp orted by thermal 
pressur e within the first four scale-heights (Kim, Hong, Ryu 
& Jones 1998). On the observational side, Rupen (1991) 



presents evidence that the non-thermal pressure gradients 
equal 400% of the thermal pressure gradient in NGC 891, 
and 50% in NGC 4565. Since the luminous mass distribu- 
tions in these galaxies are very similar, one would expect 
very different gas layer widths if non-thermal effects were 
indeed important. In fact, these two galaxies have almost 
indistinguishable flaring curves, and we can conclude that 



thermal pressure gradients dominate (Rupen 
1996aQ 

To summarize, simplified theoretical considerations in- 
dicate that non-thermal pressure support could dominate 
the hydrostatic balance if extreme assumptions about the 
magnetic field geometry and the cosmic ray energy density 
are made. More realistic assumptions regarding the vertical 
gradients of Pb and Pc, as well as observational data lead 
to the opposite conclusion. 

In the next paragraph we will investigate the effects 
of significant non-thermal pressure support by treating the 
gaseous velocity as an unknown, notwithstanding the indica- 
tions that non-thermal pressure support is actually small in 
the outer Galaxy. Based on the arguments presented above, 
we expect a' g to lie between ~0.85 and v2 times the default 
value of 9.2 kms" , while a more realistic upper limit to 
a' g /a g might be 1.05 [cf. eqn. (§§)]. 



4.4 Constraints on the pressure term 

In this section we address the question as to how the in- 
terrelations between the Galactic constants, E* and q are 
affected when assuming that a g is unknown. We follow the 
procedure to determine the halo flattening as outlined in 
the Introduction (and Paper I) for several trial values of a g 
{a g = 8.0,8.6,9.2,9.9,10.5,12 and 14 kms" 1 ), and on the 
same Ro — ©o grid as in Paper I. As predicted by equa- 
tion (|l^), the inferred halo flattening and stellar column- 
density are rather sensitive to the adopted value of a' g . In 
figure ^ we present q (dotted lines) and E* (long dashed 
lines) as calculated for increasing values of a' g , from left to 
right, and from top to bottom. The oblate region of param- 
eter space, below the heavy dash-dotted line, is strongly re- 
stricted if u'g/o-g < 1 (lower-left panel), while essentially the 
whole range of the Galactic constants is allowed if the effec- 
tive dispersion is as large as 12 kms" 1 (upper-right panel). 
From these figures we can also infer that the mass of the 
stellar disk, if measured accurately, constrains the allowed 
range of a' g . 

Values for a'„ that are as much as 30% larger than the default 
value (upper right panel of Fig. ^) are completely excluded 
because of the very small stellar disk masses required. On 
the other hand, lower velocity dispersions need lower density 
(rounder) halos and hence a more massive stellar disk, for 
given Galactic constants. This is indeed what is observed in 
the lower-left panel of Fig. ^|for Oo > 175 kms" 1 . However, 
below ©o~175 kms" , the situation is much more chaotic 
(extreme negative values for E, and q), which is the result 
of the fact that the H I flaring and the stellar kinematical 
constraints are mutually exclusive (Appendix |b|) . In fact, an 
effective velocity dispersion < 8 kms" 1 is essentially ruled 
out. 



4.5 Other Constraints 

In order to take all constraints properly into account, we 
1991' Oiling should present a three dimensional plot with Ro, ©o and a' g 



as the axes. Because this is a little tedious, we opt for to 
determine the halo flattening and E* along lines of constant 
Qo ■ We select Qo values w hich bracket the u ncertainty of the 
proper motion of SgrA * ( |Reid et al. 199S| ). In each of the 
three panels of Figure | (fi = 29.9, 27.4, 24.9 kms" 1 kpc" 1 , 
from top to bottom) we present contours of constant halo 
flattening by dotted lines as a function of Ro and a' g . The 
heavy dashed-dotted line is the round-halo contour. The 
hashed parts of the diagram depict regions of parameter 
space where E*= 35 ±5 M0 pc" 2 (heavy horizontal hash) 
and E*= 27.8 ±3.8 M Q pc" 2 (light vertical hash). These 
columndensity ranges correspond to the stellar disk mass as 
determined by KG89 and GBF97, respectively. In figure ^ we 
also plot our estimate to the upper limit of the non-thermal 
pressure support (a' g ~ 1.05a 9 ) as the thick horizontal line. 
The cross at (Ro,<y' g ) = (7.2,9.2) corresponds to the value 



of Ro derived from the Oort constant ( Oiling & Merrifield 



1998c ) and the standard value of the gaseous velocity dis- 
persion. 
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R [kpc] R [kpc] 

Figure 2. Here wc present the mutually consistent set of Galactic constants, stellar columndensity in the Solar neighborhood (£*; long 
dashed lines) and halo flattening (q; dotted lines). The individual panels show the results of the calculation for several values of o-' g /a g : 
0.85 (a' g =8 kms -1 , lower left), 1.0 (a' g =9.2 kms -1 , lower right), 1.15 (a-g=10.5 kms -1 , upper left), and 1.30 (a' g =12 kms -1 , upper 
right). The oblate-prolate boundary is indicated by the heavy dash-dotted line. The heavy full line and the heavy dashed line corresponds 
to KG89's determination of £*, and the ±1 — <r values. The upper limit of GBF97's determination of the stellar columndensity corresponds 
to the 30 Mq pc -2 contour. Both the halo flattening and stellar column arc determined to ~ 6% accuracy (see also Appendix |c[) . 

Because the Galactic constants, £*, and q have to be mutually consistent one cannot arbitrarily choose the four parameters. Fixing one 
of the four parameters severely restricts the other three. Any two parameters follow immediately from any choice of the other two, for 
a given a g . In the shaded region of parameter space, the mass of the stellar disk is as measured by KG89 (dark shading) and GBF97 
(light shading). 



A general feature of these diagrams is that flatter ha- 
los are found in regions with large values of u' g . This arises 
naturally from the fact that a stronger gravitational pull 
is needed to constrain a gas layer with additional pressure 
support, for a given observed thickness [see also eqn. jis])]. 



These figures also clearly show that regions with large DM 
densities, due to either large Qo or small q, have low stel- 
lar columndensity in the Solar neighbourhood. We can seen 
that significant non-thermal pressure support, which we de- 
fine here as having a' g — v2cr s ~ 13 kms -1 , requires very 
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6.5 V 7.5 8 8.5 

R o = o /29.9 [kpc] 




6.5 7 7.5 8 8.5 

R o = o /27.4 [kpc] 



©o 170 180 190 200 210 




R = © /24.9 [kpc] 



Figure 3. In this figure we present the interrelations between the adopted value of the gaseous velocity dispersion (erg') and the Galactic 
constants. The line coding is the same as in Figs. ^ This figure was generated by extracting the halo flattening and S* values along 
the lines @q/Rq = 24.9, 27.4, and 29.9 (from top to bottom), for several values of the effective gase ous velocity dispersion. Th e fat cross 
represents our best estimate for the Galactic const ants \Rn = 7,1 ± 0.4 kpc, ©o = 184 ± 8 kms , (Oiling & Merrifield 1998c)], and the 
gaseous velocity dispersion [<r g = 9.2 ± 1 kms -1 , (Malhotra 1995)]. The thick horizontal line represents our theoretical expectation as 
to the maximum value of <j' g (eqn. ^). Large effective velocity dispersions require flattened halos so that the model gas layer widths are 
as thin as observed. Round, and even prolate, halos are found for small a g ' values. In the shaded region of parameter space, the mass of 
the stellar disk is as measured by KG89 (horizontal dark shading) and GBF97 (vertical light shading). 
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highly flattened dark matter halos and very small stellar 
columndensities in the Solar neighbourhood. Such strong 
non-thermal pressure support is thus ruled out. 

Several other generic conclusions can be drawn from 
figure ^. First, highly flattened halos are only possible for 
realistic values of E* if _Ro < 7 kpc, and only for small val- 
ues of Qo- Second, large values for Ro, an acceptable stellar 
columndensity, and an oblate halo occur only if the angular 
velocity of the Milky Way is small, whatever the value of 
the effective velocity dispersion. Third, the extra degree of 
freedom associated with a'„ does not greatly influence the in- 
ferred halo flattening if a good constraint on E* can be used. 
And finally, the combined constraints set by the observed 
stellar columndensity and the halo's oblateness severely re- 
strict the allowed range for cr' in particular for Ro > 7 kpc. 

More specific results follow if we are willing to make 
more restrictive assumptions. For example, if we assume that 
Ro > 7 and a' g < 10.2, we find q £ 1.0 (fi = 29.9), q £ 
0.65 (fi = 27.4), and q > 0.2 (fi =24.9). In case Rq > 8 
kpc, q > 2, q > 1.3 and q > 0.8 for the same angular veloc- 
ities. Alternatively, when choosing particular values for the 
Galactic constants, simple relations between the remaining 
three parameters of the mass model follow. For example, tak- 
ing the IAU-recommended values for Ro and <do (8.5 kpc, 
and 220 kms" 1 ) we find q ~ 2.26-0.90 x da' g +0.11 x da' g 2 
and E* ~ 65.6 - 35.2 x do' g + 4.8 x da' g 2 , with da' g =a' g - 
9.2. Also, with the standard value for the gaseous veloc- 
ity dispersion we find: Oo(g < 1) <J 187 + 5 x dRo and 
E„(0.5 < q < 1) ~ 37.5 + 18.4 x dR + 8.5 x dR%, with 
dR =Ro-7.5 (cf. Fig. g the lower-right panel). 

An inspection of the lower two panels of figure ^ re- 
veals that rather tight constraints can be placed on the halo 
shape, the local angular velocity and the effective velocity 
dispersion if the halo is oblate and 30 < E* < 40 Mq pc -2 
and Ro > 7 kpc. In that case we find: 24.9 <; Qq <. 27 .4 
kms" 1 kpc" 1 , 0.5 < q < 1, and 8.6 < a' g < 10.3 kms -1 . 
If we impose the additional constraint that the non-thermal 
pressure support is limited to 5% as derived in equation (p5|), 
it follows that the Sun's distance to the Galactic centre is 
less than 8 kpc, and that the rotation speed of the Milky 
Way is less than 200 kms -1 at the Solar circle. 

The stellar disk mass provides a strong constraint on 
the effective velocity dispersion of the gas: a low disk mass 
requires more dark matter, which would lead to a thinner 
gas layer in the outer Galaxy if a' g were not increased. For 
example, if the stellar disk mass exceeds 22 Mq pc" 2 , fig- 
ure ^ shows that the effective velocity dispersion has an up- 
per bound of about 10.5 kms -1 , so that the non-thermal 
pressure support can not exceed 14%. 



5 SUMMARY AND CONCLUSIONS 

In Paper I we showed that the constraints on the dark matter 
density in the Solar neighbourhood and the outer Galaxy 
place tight limits on the choice of the Galactic constants, 
the mass of the stellar disk, and the dark halo's flattening. 
The internal errors for this procedure are of order 6% for 



both the halo flattening and the local stellar column density 
(see Appendix ^). 

In this paper we present the analytical tools that al- 
low for an efficient search through parameter space. Em- 
ploying these tools, we extend the analysis of Paper I by 
investigating the effects of the ill-determined contribution 
of non-thermal pressure support on the H I flaring, and the 
consequences for the inferred Rq, Qq, S*, and q values. The 
strongest constraints available are the observed Oo /Ro ratio, 
the mass of the stellar disk and the fact that the dark halo 
is almost certainly oblate. Taken together, these constraints 
rule out substantial contributions to the support of the H I 
layer by cosmic ray pressure or magnetic fields if the Sun's 
distance to the Galactic centre is greater or equal than 7 
kpc, in agreement with theoretical predictions. We find that 
the Milky Way's dark matter halo is close to spherical for 
all but the smallest values of _Ro or Go. A dark matter halo 
as flattened as q = 0.2 is only possible if our distance to the 
Galactic centre is smaller than about 6.8 kpc. 

It is possible to construct a self-consistent oblate model 
of the Galaxy with 7?o = 8.5 kpc and Oo = 220 kms" 1 , but 
only if the local stellar columndensity is less than about 18 
Mq pc -2 , and a' g J> 11 km s _1 . 

Kuijken & Gilmore's (1991) determination of the colum- 
ndensity of matter within 1.1 kpc of the plane (71 ±6 
Mq pc -2 ) is robust and valid over a wide range of Galactic 
constants and disk scale- lengths. For the stellar contribution 
to this total mass we suggest a consensus average of E* = 35 
and a consensus error of 10 Mq pc -2 . 

If Ro > 7 kpc, then the dark halo of the Milky Way is 
fairly close to spherical, independent of the amount of non- 
thermal pressure support. Such a round halo argues against 
dissipational baryons as a viable dark matter candidate. Fur- 
ther, since all other baryonic dark mat ter candidates have a l- 
ready been observationally excluded ( Hegyi fe Olive 1986 ), 



we must conclude that the dark matter halo of the Milky 
Way is most likely made up of something altogether more 
exotic (MACHOs, neutrinos, axions, neutralinos ...). Even 
in the Solar neighborhood, there are non-negligible quanti- 
ties of this material: our proposed dark halo models imply 
that it amounts to some 0.42 GeV/c 2 per cubic centimetre 
or (11 ±5) oiMq pc -3 . The direct detection of this material 
remains a challenge for experimental physicists. 

Employing the local disk mass and the flaring of the 
Galactic H I layer, we find strong correlations between the 
parameters of axisymmetric mass models. We presented sev- 
eral examples in the previous section. Since we make strong 
predictions as to the values of Ro, Oo, E„, q and a' g (Figs. PJ 
and ^|), our models can be subjected to experimental verifi- 
cation. For example, all parameters but a' g could be deter- 
mined using astrometric data from future astrometric space 
missions such as FAME, SIM and GAIA. Such high precision 
data are ideally suited to support, or falsify, the models we 
propose here. It would also be worthwhile to investigate the 
effects of deviations from axisymmetry on the inferred halo 
shape and rotation speed. Whatever the outcome, we will 
learn a great deal more about the structure and dynamics 
of the Milky Way galaxy. 
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APPENDIX A: MORE ON THE VERTICAL 
DISK-HALO CONSPIRACY 



In section 3.1 we reviewed the stellar kinematical route to 



determine the local stellar columndensity. In this appendix 
we try to emulate and check this method in some detail. 
In order to do so we calculate the exact vertical force law 
-Kz, exact for a limited number of Galaxy models listed in 
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Table ^| using the global approach (§ |3,3| cf. Oiling 1995). 
We also compute an app roximation to Kz^act from the local 
mass distribution (§ 3.3; cf. Oiling 1995), which follows from 



an integration of the Poisson equation: 



(z) = -AnG I dz [p m (z') + Prot 



with p m the matter density, p ro t{= 



(Al) 



a pseudo 



AttG R dR 

density which arises from gradients in the rotation curve 
and K z the v erti cal force per unit mass 
In figure 



Al 



we present true vertical force (filled circles) 
and the local approximation (open circles) for one particular 
combination of parameters which lie close to the best esti- 
mates for the Milky Way (|Kuijken fc Gilmore 1989, 1991 



Oiling & Merrifield 1998c). We see that in the Solar neigh- 
borhood and close to the plane, the local approximation is 

reasonably accurate^] and can be used as an approximation 
to the true vertical force law. Significant deviations from the 
true K z only occur above approximately 800 pc. We also in- 
clude the observationally determined K z [crosses with error 
bars; KG89's equation (1) with D=250 pc]. 

Now we investigate the possibility that the actual dark 
matter density law differs from the distribution in our mod- 
els [eqn. (jjj)]. Choosing a different DM distribution while 
keeping the observed K z fixed requires a different lumi- 
nous mass distribution. As a simple test case, we plot a 
model where the DM de nsit y is multiplied by 1.92 (labeled 
"dGR+Hxl.92" in Fig. Al). In this case, we have to de- 



crease the stellar columndensity to the GBF97 value (25.8 
Mq pc -2 ) in order to keep K z approximately unchanged. 
Next we consider DM distributions inspired by Bahcall's 
"P" models (1984b), where the DM density is assumed to 
be a linear combination of the known components: 



p err = Xd X pd + X g X p g + Xh X ph + Xrot X p r ot , (A2) 

where the subscripts d, g, h denote the disk, gas, and halo 
components, respectively. Using this formalism, we can cal- 
culate the columndensities of all components (ISM, stars, 
DM and Ej^) which arise from a particular choice of the Xi- 
values. From such an erroneous density distributions (p err ) 
we deter mine the erroneous vertical force (K z ^ rr ) using 
equation (Al). We present the K z .true/ ' K z , err ratio, aver- 
aged in several z-height ranges, in Table Al (columns 5-7). 
In this table we also list the Xi, E* and E tot values for some 
models. In addition to the P-models, we also evaluate models 
in which we keep the ISM and p ro t contributions fixed while 
varying the stellar and dark matter densities such that the 
force at 400 pc above the plane equals the true value. From 



II For the ensemble of models for which the parameters differ by 
< 1-cr from t he m odel presented here, the approximations based 
on equation ( |Al| ) reproduce the true K z typically to within 7%: 
sometimes models that include the rotation curve gradient term 
perform best, sometimes not. The same is true for the models 
listed in Table U with other Galactic constants than those pre- 
sented in figure] 



Table Al. The accuracy to which various erroneous density dis- 
tributions can reproduce the true vertical force. Th e first four 
columns describe which p(z) was used [cf. eqn. (A2)]. The next 
three columns indicate the average ratio of Kz,true/^z. err Over 
the vertical ranges indicated [0 - 0.4 kpc (Zl), 0.4 - 0.8 kpc (Z2), 
and 0.8 - 1.2 kpc (Z3)]. The typical variation within these radial 
ranges is a few %, except for the gas-only case where it can be 
up to 35%. We also indicate the model values for S* and ^\ a \, in 
units of Mq pc - 2 . The last two lines give the average and stan- 
dard deviations of the lines above (excluding the gas-only model) . 



x d 




X ro t 


x h 


Zl 


Z2 


Z3 


£* 


yl.l 
U tot 


1.25 


1.25 


1.25 


0.00 


1.03 


0.99 


0.93 


45.4 


63.6 


1.29 


1.29 


0.00 


0.00 


0.95 


1.04 


1.17 


46.9 


65.6 


2.05 


0.00 


0.00 


0.00 


0.85 


1.03 


1.01 


74.5 


74.5 


1.95 


0.00 


1.95 


0.00 


0.85 


1.03 


1.04 


70.8 


70.8 


0.00 


3.49 


0.00 


0.00 


1.35 


0.90 


0.72 


0.0 


50.6 


1.40 


1.00 


1.00 


0.03 


1.03 


0.99 


0.93 


50.9 


66.1 


1.25 


1.00 


1.00 


0.46 


1.03 


0.99 


0.93 


45.4 


70.6 


1.00 


1.00 


1.00 


1.17 


0.97 


1.02 


1.08 


36.3 


78.1 


0.74 


1.00 


1.00 


1.92 


0.95 


1.04 


1.17 


26.9 


86.1 


0.50 


1.00 


1.00 


2.59 


0.93 


1.06 


1.24 


18.2 


92.9 


0.25 


1.00 


1.00 


3.30 


0.92 


1.08 


1.33 


9.1 


100.4 


0.12 


1.00 


1.00 


3.67 


0.91 


1.09 


1.37 


4.4 


104.3 



39± 

58% 



77± 
21% 



table Al we see that it is possible to decrease E* by more 
than an order of magnitude and simultaneously increase Xh 
without altering -?Cz(|z| < 800pc) substantially. It seems that 
the vertical distribution of dark and stellar matter conspire 
in such a way that their relative contributions can not be 
easily separated. Thus, this "vertical disk-halo conspiracy" 
is much like the classical disk-halo conspiracy which arises 
in considerations of the radial distribution of luminous and 
dark luminous matter derived from gal axy rotation curves. 

it is clear that, 



Al 



From the data presented in Table 
without additional assumptions, stellar kinematical data ex- 
tending to a few hundred parsec are not sufficient to de- 
termine the mass of the stellar disk. However, Y,\ n \ is much 
better determined (rightmost column of Table |Al| ), in agree- 
ment with Kuijken & Gilmore (1991). In order to discrimi- 
nate between models, it is essential to incorporate a rotation 
curve constraint to limit the possible DM densities. Kuijken 
& Gilmore (1989) used such a constraint explicitly, while 
Bahcall (1984b) and Flynn & Fuchs (1994) used "reason- 
able" values for the density of non-disk-like dark matter to 
rule out extreme values for the mass of the stellar disk. Fur- 
thermore, if data at larger z are consid ered , the differences 
between the force laws shown in figure Al start to become 
apparent and can be used to limit the allowed range for E* 
(Flynn & Fuchs 1994). 

Although th e th e model parameters we have used to 
construct Figure Al and Table Al lie at the extreme end 
of the models investigated by Kuijken & Gilmore, we find 
that other models, within the l-2cr range of Qo/Ro , ^-d, E*, 
and Ej i, yield similar results. Thus, mass models for the 
Milky Way have to conform to the constraint T,j \— (71 ±6) 
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Figure Al. We present the exact vertical force (filled circles) at Galactocentric distance Ro for a model that almost reproduces the 
values for E* and Sj Q J proposed by Kuijken & Gilmore (1989, 1991). The local approximation to the exact force law that uses the slope of 
the rotation curve (open circles) or assumes a flat rotation curve (dashed open circles) are also indicated. Some other models in which the 
DM is distributed like various superpositions of the known components are also plotted: 1) puM = 1-29 x (p* + p g ) (labeled "DGnxl.29'); 
2) Pdm = 1-25 X (p, +p g +Prot) ("DGRxl.25"); 3) p DM = 2.05 X p, ("Dnnx2.05"); 4) p DM = 1.95 X (p, + p rot ) l"'DnRxl. 95"1: 5) 
Pdm = 3.49 X p f"nG nx3.49"). For the curve labeled "dGR+Hxl.92" we used a smaller columndensity of stars 25.8 Mq pc~ 2 (Gould 
Bahcall fe Flynn 1997), an d mu ltiplied the halo contribution by 1.92. Most of the models have very similar force laws in the region close 
to the plane (see also Table [Al] ). We also include KG89's vertical force law and the errors thereupon (crosses with error bars). 



M pc _2 .To summarize the above, in order to determine the 
stellar columndensity from stellar kinematics, the existence 
of the vertical disk-halo conspiracy requires that one has 
to include self-consistent rotation curve constraints and/or 
sample the region above 800 pc. 



APPENDIX B: STELLAR KINEMATICS AND 
THE RELATION BETWEEN THE MODEL'S 
PARAMETERS 

In section Vi/A we described how the observationally de- 
termined E tot -value imposes correlations between luminous 
and dark matter in the Solar neighbourhood. In this sec- 
tion we present some specific examples to illustrate what we 
can learn about the structure of the Milky Way by treat- 
ing E* as a parameter which is only constrained by E^J [cf. 
eqn Note that these correlations are independent of any 
H I flaring constraints. 



At the Solar circle, and for z — 1.1 kpc, equation 
can be rearranged to read: 



;l.LRo +Ro (0.9-R )q 



-13/20 



q<0.15 (Bl) 
g>0.15 (B2) 



where the core radius has been set to Ro- In this appendix 
we use -for the sake of convenience- the symbol q for the 
halo flattening derived from the local stellar kinematics con- 
straint, and gjr j for the halo's shape derived from the H I 
flaring. Taking the fiducial values for Ej J, E», and E g to 
be 71, 35, and 14.5 Mq pc -2 , respectively, and employing 
equation (^) we determine the fiducial value for E^' 1 to be 
21.5 Af pc" 2 . Further, defining E, = 35 + <5E«, use the 
fiducial columndensities defined above, and apply equations 
(Bl) and (B2), we can relate the required halo flattening to 



<5E„ 



21.5 -Ro E^fl) - 3S» 
0ARpF(l) 



(B3) 
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yl.l 



(1) 



21.5 



20/13 



1 



/ 20/13 
I 21.5 



<5E» 



(B4) 



for the same q-ranges as in Eqns. (Bl) and (|B2j). Since 
Ej l ' 1 (g = 1) depends on 0o [eqn. (Q)], the solutions of equa- 
tions (B3) and (B4) depend on Rq, 0q, and <5E«. We use our 
mass models to determine Hi' 1 (a = l;i?o,0o) and rewrite 



<5E* = we find: 



equations (B3) and (B4) for a few interesting cases. For 



q p 


s 0.07 ±0.05 


(8 


= 165) 


(B5) 


q p 


d 0.4-0.1(Rq -7.1) ±0.1 


(00 


= 175) 


(B6) 


q p 


s 0.8-0.2(i? -7.1) ±0.2 


(0o 


= 185) 


(B7) 



Thus, lower rotation speeds and larger -Ro's require more 
highly flattened halos. Furthermore, for any fixed value of 
Ro, the slope of the q — <5E* relation depends strongly upon 
0o. For example, with Tio— 7.1 kpc we find: 



q ^ 0.069 + 0.005 <5E» 
q ^ 0.400 + 0.030 (5E„ 
q ^ 0.800 + 0.060 5E* 



(0o = 165) 
(0o = 175) 
(0o = 185) 



(B8) 
(B9) 
(B10) 



Equations (B5) through (B10) clearly show that small val- 



ues of the Galactic rotation speed imply a highly flattened 
dark matter halo, whatever our distance to the Galactic cen- 
tre, and whatever the mass of the stellar disk. We also see 
that the last relations constrain E* rather tightly: for mod- 
els with 0o ~ 185 (175) kms , an increase in the stellar 
columndensity of ~ 3 (20) M@ pc~ 2 covers the whole allowed 
range for q. 



In section 4.4 and figure y we have see the stellar kine- 
matics and H I flaring constraints provide mutually exclu- 
sive constraints on E, and q if both the rotation speed and 
the effective velocity dispersion are s mall. This can be un- 
derstood as follows: equations (p35|)-(BlC) show that the 
E tot constraint implies small q- values for low 0o's. Like- 
wise, the H I flaring constraint yields small q's for small 
Go's, but only if the velocity dispersion does not decrease 
by too much [cf. eqn. (|l6|)]. Furthermore, because the slope 
of the g(E«) relation becomes shallower with decreasing 0o, 
the accessible range for gi.i decreases with 0q. If the flaring 
analysis leads to a gjj j ~ 1 halo, extreme values for E» are 
required to match qn with gjj j . Obviously, if the required 
stellar column exceeds T,j \, it is not possible to construct 
a self-consistent model for that particular combination of 
Galactic constants and gaseous velocity dispersion. In these 
circumstances it is possible that our procedure to determine 
E» and q yields extreme values for E* and q, and sometimes 
even negative values. 

APPENDIX C: ERROR ESTIMATION 

At this point it is also possible to estimate the accuracy to 
which the various parameters in figure ^| are determined. For 
example, as we have seen before, the halo flattening inferred 
from the H I flaring has only a slight dependency on the 
mass of the stellar disk (cf. figure 3 of Paper I) . This allows 
us to estimate the flattening of the halo by averaging the 



gjj t values from the various model runs at a given _Ro and 
0o. For example, the nine <?jt j errors tabulated in table ^ 
yield Sq^ j /qjr I ~ 0.06 for the two low Rq values at q ~ 0.7 
Employing q oc Oq [cf. eqn. (DjT)] we find: 



50o 
0o 



s m i 



(Cl) 



2 <7H I 

The qn j values as determined from the flaring measure- 
ments are thus precise to about 6%, while the value of the 
local Galactic rotation speed we derive from gjj j measure- 
ment has an estimated accuracy of 3 percent. 

To estimate how well the local stellar column density is 
determined from the flaring and the observed E^;i value of 



the total column density we re- write equation (B4) to read 
q = ci(l + C2<5E„), and find: 



g 

ClC 2 



(^) 2 + (-) 



5ci . 

Cl 



(C2) 



L (l)/21.5, 



C2 



(20/13)/21.5 and q/{c x c 2 ) 



0.8 and E^' (1)=21.5. If assume that ci is with- 
vanishes, and we arrive 
of about 0.7 M P "-~ 2 



5(<SE„) r 

where ci = 
11.2 for q 

out error, the second term in eqn. CI2 
at a lower limit for the error in E* ot about U.Y Mg pc 
If we assign the full error in the observed total column den- 
sity (6 M Q pc" 2 ) to E/ [ 1 (l), we have a ~ 1 ± 0.43, so 
that <5(5E*)) = 11.2^(0.06)2 + (0.43) 2 ~ 4.8M Q p C - 2 . We 
thus estimate that our method of deriving the stellar col- 
umn density has an accuracy somewhere between 0.7 and 
4.8 Mq pc -2 or about to 2 to 14 percent. 

The error estimates above are close to the errors derived 
from our detailed modeling procedure. 



APPENDIX D: 
TABULATED 



THE ISM OF THE MILKY WAY 



In this appendix we present a tabulated version of the radial 
variation of the atomic and molecular hydrogen, as well as 
their widths. We present these data in a manner which is 
independent of the values of the Galactic constants as well 
as the shape of the rotation curve. The H I columndensity at 
the Solar position includes 1.4 Mq pc -2 of ionized hydrogen. 
No ionized hydrogen is included at any other radius. The 
columndensities listed do not include the contribution due 
to Helium. In our model calculations of the potential we 
increase the listed columndensities to include 23.8% Helium. 
We have brought the data from the sources listed in section 



2.1 onto a common distance scale defined by Merrifield's 
(1992) determination of the W(R/Ro) = v ra d/(sin£cosb) 
curve. Here v ra d and £ and b are the radial velocity and 
the Galactic coordinates, respectively. In practice this re- 
scaling works as follows: 1) from the rotation curve [O'(ii')] 
in the original reference, determine R' /R' and W' (R' / R' ) = 
R' [Q'(R')/R' - O'o/R'o] [the primed quantities refer to the 
values assumed in the reference], 2) for each property X' 
(gas layer width and volume density), determine X'(W') 
from X'(R') and W'(R'), 3) find Merrifield's R/R values 
for which W = W, 4) the property X' re-gridded onto 
Merrifield's distance scale is now given by X(R/Ro). 
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Table Dl. The surface density and thickness of the atomic and 
molecular hydrogen components in the Milky Way as a function 
of scaled Galactoccntric radius, R/Rq . The H I and H2 colum- 
ndensities have units Mq pc~ 2 . The widths listed are full widths 
at half maximum (FWHM) and have units of i?o- For R > Rq, 
the widths of the H I layer is the average of several works, see 
Paper I for details. Some entries have been left blank: we were un- 
able to derive reliable widths at these radii. The columndcnsitics 
at these radii are from Woutcrloot et al. (1990). 



R/Ro 


E H I 


E H 2 


%I 


6W Rl 


w n 2 


0.218 


0.59 


1.08 


0.0330 


0.0014 


0.0151 


0.276 


1.04 


1.54 


0.0320 


0.0031 


0.0158 


0.309 


1.23 


1.47 


0.0323 


0.0017 


0.0134 


0.343 


1.30 


1.67 


0.0306 


0.0014 


0.0119 


0.376 


1.29 


2.67 


0.0294 


0.0008 


0.0131 


0.409 


1.31 


4.29 


0.0268 


0.0045 


0.0140 


0.440 


1.41 


4.91 


0.0290 


0.0031 


0.0142 


0.467 


1.59 


5.02 


0.0313 


0.0028 


0.0143 


0.499 


1.85 


5.57 


0.0278 


0.0011 


0.0153 


0.530 


2.10 


6.23 


0.0240 


0.0017 


0.0165 


0.559 


2.30 


5.87 


0.0301 


0.0014 


0.0168 


0.586 


2.48 


5.48 


0.0346 


0.0014 


0.0163 


0.612 


2.90 


5.18 


0.0344 


0.0023 


0.0152 


0.641 


3.72 


4.62 


0.0415 


0.0034 


0.0139 


0.668 


4.36 


3.94 


0.0504 


0.0017 


0.0129 


0.692 


4.24 


3.73 


0.0525 


0.0014 


0.0120 


0.717 


3.60 


3.97 


0.0447 


0.0031 


0.0115 


0.742 


3.17 


4.27 


0.0452 


0.0039 


0.0122 


0.764 


3.23 


4.33 


0.0485 


0.0130 


0.0136 


0.786 


3.68 


4.09 


0.0495 


0.0101 


0.0152 


0.806 


4.22 


3.60 


0.0565 


0.0039 


0.0161 


0.828 


4.72 


2.95 


0.0530 


0.0048 


0.0159 


0.848 


4.90 


2.44 


0.0570 


0.0034 


0.0144 


0.865 


4.80 


2.06 


0.0572 


0.0011 


0.0134 


0.882 


4.52 


1.70 


0.0509 


0.0014 


0.0136 


0.897 


4.41 


1.45 


0.0525 


0.0020 


0.0146 


0.913 


4.75 


1.30 


0.0523 


0.0036 


0.0158 


0.926 


5.39 


1.26 


0.0652 


0.0017 


0.0167 


0.938 


6.13 


1.29 


0.0730 


0.0020 


0.0174 


1.000 


9.25 


1.80 


0.0577 


0.0042 


0.0199 


1.058 


8.57 


1.80 








1.089 


7.77 


1.28 


0.0611 


0.0044 


0.0238 


1.193 


6.09 


0.44 








1.332 


6.56 


1.62 


0.0734 


0.0061 


0.0385 


1.457 


7.45 


1.47 








1.586 


7.62 


1.10 


0.0801 


0.0062 


0.0458 


1.721 


6.02 


0.54 








1.848 


5.30 


0.37 


0.1108 


0.0099 


0.0742 


1.995 


4.36 


0.20 








2.101 


3.41 


0.11 


0.1391 


0.0115 


0.0767 


2.177 


3.31 


0.14 








2.262 


2.84 


0.12 








2.355 


2.36 


0.09 


0.1666 


0.0143 




2.443 


2.01 


0.05 








2.532 


1.66 


0.01 
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